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Abstract 

We apply the Liouville-von Neumann (LvN) approach to open systems to 
describe the nonequilibrium quantum evolution. The Liouville-von Neumann 
approach is a unified method that can be applied to both time-independent 
(closed) and time-dependent (open) systems and to both equilibrium and 
nonequilibrium systems. We study the nonequilibrium quantum evolution of 
oscillator models for open boson and fermion systems 



Key Words: Liouville-von Neumann approach, invariant operator, non- 
equilibrium quantum evolution, open boson system, open fermion system 

PACS: 03.65.-w, 05.30.-d 



Typeset using REVTeX 



* Contributed to the journal "Condensed Matter Physics". 



1 



I. INTRODUCTION 



The system of an open system interacts with an environment and depends on time im- 
plicitly or explicitly. A closed system can depend on time only when it is placed in an 
expanding Universe. It is the parametric interaction that forces the coupling constants of 
the system to depend on time. The time-dependent system can not maintain the equilibrium 
state since the density operator defined in terms of the Hamiltonian does not satisfies the 
Liouville-von Neumann equation. Hence, the open system requires the nonequilibrium evo- 
lution through the interaction with the environment. Recently the nonequilibrium evolution 
of open quantum systems has been an issue of much concern and debate [p]]. 

There have been two widely used methods to treat the nonequilibrium quantum sys- 
tems. One is time-path integral method where path integrals are carried out over a contour 
of complex (imaginary) time to take into account the thermal nature of the system (for 
review and references, see [0]). The other is thermal field dynamics where the system of the 
open system interacts and makes a closed system with a fictitious system (for review and 
references, see ||). 

It is the purpose of this paper to apply a recently developed field theoretical method, the 
so-called the Liouville-von Neumann (LvN) approach to nonequilibrium quantum systems 
(for references, see [|J). It is a canonical method where the evolution of open systems is 
governed by the Schwinger-Tomonaga equation or the functional Schrodinger equation. The 
LvN approach is well suited to describe the nonequilibrium quantum processes because the 
same LvN equation is used both to construct the Hilbert space of quantum states and to find 
the density operators. Also the method is nonperturbative in that the lowest approximation 
is the time-dependent Gaussian wave functional for a self-interacting system. 

The organization of this paper is as follows. In Sec. II we review the LvN approach and 
apply to a time- dependent oscillator. In Sec. Ill we study three oscillator models for open 
boson systems and find their nonequilibrium evolution. In Sec. IV we study open fermion 
systems motivated by a Dirac field interacting with an external source and a scalar field. In 
Sec. V we conclude the paper. 



II. LIOUVILLE-VON NEUMANN (LVN) APPROACH 

The LvN approach is based on the idea by Lewis and Riesenfeld that the LvN equation 
can be used to find the exact quantum states of a time-dependent quantum system M. 
In particular, they solved exactly the time-dependent quantum oscillator in terms of an 
auxiliary equation. Recently, the author has found an oscillator representation, where the 
Fock space is constructed in terms of the classical solution HJ6|-|9|]. Furthermore, the selection 
rule was put forth for the vacuum state for the time-dependent oscillator, which satisfies all 
criteria frequently imposed for well-known systems [H. 

We elucidate the LvN approach to the nonequilibrium quantum evolution of open quan- 
tum systems. The LvN approach is based on two assumptions which are physically well- 
grounded. The one is that an open system should be described by a time- dependent quantum 
system obeying the Schrodinger equation (K — 1) 

4l*(*)> = #(*)!*(*)>■ C 1 ) 
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The other is that the statistical nature of an ensemble of the system should be described 
by the density operator satisfying the quantum LvN equation. In methodology the LvN 
approach is based on the the following two theorems. 

Theorem I. If 0(t) is a time-dependent invariant operator satisfying the quantum LvN 
equation 



ifd(t)=i^6(t) + [6(t),H(t)] = o, 

then its eigenstates 

6{t)\X n ,t) = X n \X n ,t) 
provide the exact quantum state for Eq. ([[]) by 

W)> =E c « ex p[ 2 / dt(X n ,t\i—- H{t)\X n ,t)]\X n ,t). 



(2) 



(3) 



(4) 



This theorem was first proved by Lewis and Riesenfeld. In order to prove the theorem we 
use Eqs. ([D and (§) to get the identity 



(6{t) - xj 



i^-H(t))\X n ,t) 



0. 



Then Eq. ([5]) implies that 







(i—-H(tj)\X n ,t) =a n (t)\X n ,t). 



(5) 



(6) 



Theorem II. If Oi(t) and 2 (t) satisfy the LvN equation 
the same LvN equation. This follows from the identity 



then 6 1 (t)d 2 (t) satisfies 







(dxi^it)) + [(6x^)6^)), H(t) 



{iT-o^t) + [d 1 (t),H(t)]}d 2 (t) + 6 1 (t){i^-d 2 (t) + [d 2 (t),H(t)]} = o. 



o 



dt 



dt 



(7) 



In general, one can show that F[0(t)} satisfies the LvN equation, if it is an analytical 
functional of 0{t) and 0{t) satisfies the LvN equation. 

A few remarks should be made. Firstly, the LvN approach can also be applied to time- 
independent systems without any modification. In this case one can choose O = H and 
the energy eigenstates are the exact quantum states as expected from quantum mechanics. 
Secondly, according to the theorem II any invariant operator 0(t) can be used to define a 
density operator 



Po{t) 



1 



Tr(e-^W) 



-pd{t) 



(8) 



To describe the initial thermal equilibrium state a particular kind of invariant operator 
will be selected. In this sense the LvN approach can be applied to the non-equilibrium 
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quantum evolution of time-dependent systems. Finally, the LvN approach can be extended 
to quantum fields whose coupling constants depend on time. When a field is appropriately 
mode-decomposed, the system of the quantum field becomes an infinite system of finite- 
dimensional quantum mechanical systems. 

We now turn to the most typical system of time-dependent oscillators. These systems 
occur in all areas of physics and play a prominent role partly because any system around 
an extremum can be approximated by a system of oscillators and partly because the clas- 
sical and quantum evolution can be found exactly. Even a free field can be treated as a 
system of infinite number of oscillators via mode-decomposition. Thus the study of a time- 
dependent oscillator will shed light on the nonequilibrium quantum evolution of a generic 
time-dependent quantum system. 

Without loss of generality, we consider an oscillator whose mass and frequency depend 
on time explicitly: 

In Refs. m| a pair of invariant operators are found 

A(t) = i(u*(t)p - m{t)u*(t)q) , 

A*(t) = -i(u(t)p-m(t)u(t)q\ (10) 
where u satisfies the classical equation of motion 

u(t) + ^Qu(t) + u\t)u(t) = 0. (11) 

These invariant operators can be made the annihilation and creation operators by imposing 
the standard commutation relation for all times 

[A(*),i t (t)] = l, (12) 

which is guaranteed by the wronskian condition 

m{t)(u*{t)u{t) - u(t)u*(t)) = i. (13) 

From the theorem II we can find a particularly simple functional of two invariant oper- 
ators fllOf ) called the number operator 

N(t) = AHt)A(t). (14) 

The vacuum state and the number states are defined by 

A(t)\Q,t) = 0, 

|n,t) = -L(At(t))|0,i>. (15) 



According to the theorem I the exact quantum state (U) is a linear superposition of number 
states. The density operator is defined by 
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p(t) = le-^W, (16) 



where Z is the partition function 



Z =J2(n,t\e-^ oN ^\n,t) 

n=0 



= 13^' ( 17 ) 

The density operator has the coordinate representation 

1 oo 

p(x,x',t) = - (x\n,t)(n,t\e-^" A '^ A ^\n f ,t)(n f ,t\x) 

n,n'=0 
i oo 

= zT,**( x > t )KW,t)e- fluon , (18) 

n=0 

where 

*„(x,t) = <x|ra,t) (19) 
are the harmonic oscillator wave functions. 



III. OPEN BOSON SYSTEM 

In this section we shall consider three oscillator models for open boson systems. These 
systems are motivated by a scalar field in a linear or quadratic interaction with an envi- 
ronment. The first two models have all the essential characteristics of a scalar field with a 
quadratic interaction with an external source 

i Tfl^ 1 

C = -0^(x,t)0^(x,f) - -f<P 2 ^t) + j(x,t)<?(x,t). (20) 
The last model mimics the scalar field interacting linearly with the external source 

1 777^ 

C = -«9 M 0(x, t)<9^(x, t) - ^0 2 (x, t) + j(x, t)0(x, t). (21) 

As the simplest open system we consider a time-dependent oscillator described by the 
single- mode Hamiltonian 

H(t) = n D (t) (ata + \) + n<°> (t) (\a* + la 2 ) , (22) 

where 

a = Vr « - (23) 
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These operators are the standard annihilation and creation operators of the time- 
independent oscillator 

H = \p 2 + fg 2 - (24) 
We look for a pair of first order invariant operators of the form 



/ A(t) \ _ ( u{t) v(t) \fa 



(25) 



where 

l Jt{v(t)) + {nW(t) -n( fl )(t)JU(*)/ { } 

And we further impose the standard commutation relation 

[A(t),A*(t)) = l, (27) 
which leads to the unimodular matrix 



v*(t) u*{t) 
The inverse transformation is given by 



a\_(u*{t) -v{t)\(A{t)\ 



With respect to the vacuum state ( fL5f) the position and momentum dispersion relations 
are calculated to be 

(f)(vac) = 77— K0 - KOI 2 ) 

(^<t™) = y|t*(f) + w(f)^ (30) 
and the uncertainty relation to be 

(q 2 ) (vac) (p 2 ) (vac) = \{ l ~ («(*)«*(*) - u*(t)v(t)) 2 }. (31) 
By parameterizing Eq. (28|) 

= cosbW", u(i) = sinhre^", (32) 
the uncertainty relation can be rewritten as 

if) (vac) (P 2 ) (vac) = J{l + sinh 2 r sin 2 (6> u - 9 V )}. (33) 
Using the density operator (|T6|) , the thermal dispersion function is found 
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(tfW) = ^-coth(^)\u(t)-v(t)\". (34) 



The second model is a system of coupled oscillators whose coupling constants depend on 
time. The Hamiltonian is described by 

Hit) = Eng? + bap) + n%\t)hlal + rt$\t)\a a a Pl (35) 

af3 " 

where VL* a/3 = Vtp a to guarantee the hermiticity of the Hamiltonian. As for the first model 
we may find pairs of first order invariant operators 

Aa(t) = E u ap{t)ap + Vapitfajj, 
P 

4(*)=E<^(*)^+</3(*)4' ( 36 ) 

p 

It is straightforward to see that Eq. (Q) leads to the system of equations 

iu a p(t) + E d${t) Uai - EK?*(*) + rt°>(t))v a7 (t) = o, 

7 7 

iv a p{t) - E^? } (*K 7 + EK>) + ^(*))^ 7 (*) = 0. (37) 

7 7 

The commutation relations 

[A a (t),A^t)]=6 a p, (38) 

lead to the conditions 

E( u *y(*)«/W*) ~ v *y(*) u i»r(*)) = <W- (39) 

7 

The final model for the open system is a system oscillator coupled to an environment 
described by the Hamiltonian 

H = H( sys ) + H( sys - env ), (40) 
where the system and the interaction are 

H( S ys) = uj a^a, 

H {sys _ env) = u(t)a ] + v*(t)a. (41) 

Here, u(t) is a complex function and of particular interest is a stochastic function. Extending 
the algebra by including the identity operator e, a pair of first order invariant operators of 
the form are found 



A(t) \ _ ( u{t) v{t) \fa\f w{t)e 
A\t)J \v*{t) u*{t))\a)) + W{t)e)- 

The LvN equation (0) is satisfied with the solutions 



(42) 
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u(t) = u e iuJot , v(t) = v e- iuJ0t , 

w (t) = i J* dt(u v{t)e iuJot - v u*{t)e- iuJot ). (43) 



The Fock space is similarly constructed as in Eq. (|15| 
Using the inverse transformation 



u*(t) -v(t) \ f A(t) -w(t)e 
at) ~ {-v*{t) u{t) J \A\t)-w*{t)e 



(44) 

it is straightforward to obtain the thermal expectation value of the system Hamiltonian 

(H(sys))(ther) = U (u* Uo + VqVq) ^ _ ^ + ^O^O^O + U \u(t)w*(t) - v(t)w(t)\ 2 . (45) 

The stochastic force has the correlation functions 

(u(t)) M = (u*(t)) {u) = 0, 

(v*(ty(t)) M = f(t'-t), (46) 

where the average is taken with respect to a probability distribution of fluctuating parame- 
ters. For a delta-correlated force 

(if(f)v(t)) {v) =D6V-t), (47) 

the energy expectation value becomes 

(H(sys))(ther) = ^0 («Q M + VqVq) ^ _ ^ + ^0VqV + LO Dt. (48) 

Thus the expectation value increases in proportion to the time of interaction with the ex- 
ternal driving force (0). Similarly, for a Gaussian-correlated force 

(v*(t')v(t)) (v) = D ie -^-\ (49) 

we obtain the energy expectation value 



(H( sys ))(ther) = U (u* Q U + W ) ^ _ - + UqVqVq 

2Dj 2 2D-fe-^ 



1 , 2^7(72 - ^ 

-5— + UJ V V S- o 

eP Uo - 1 7 +^o 



0/ 



H — 5 otH 5- — n" 7 - u)q) cos ujQt-2^ujQSVnujQt\. (50) 

The expectation value has constant terms, a linearly increasing term and oscillating but 
decaying terms. 

We are able to get the well-known results for these models [|]]. However, the LvN ap- 
proach provides us with a clear picture for the open quantum system by making use of all the 
lessons from quantum mechanics and proves technically much simpler and straightforward 
than other methods. 
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IV. OPEN FERMION SYSTEM 



It was not until a quite recent time that the LvN approach has been applied to fermion 
systems JTD , |TT | . The description of fermion systems requires not only the annihilation and 
creation operators b, b' for the particle but also those d, dJ for the antiparticle. These oper- 
ators satisfy the standard anticommutation relations 

{S, P} = {d,$} = i, 

{S, b} = {ft, ft} = {d, d} = {S, $} = 0, 

{6, d} = {b, = {ft, d} = {ft, $} = 0. (51) 

The interaction between the fermion system and an external source requires the Grassmann 
variables c and that satisfy the anticommutation relations 

{c,ft} = l,{c, £} = {c t ,c t } = 0, 
{c,b} = {c,ft} = {c,$ = {cj} = 0, 

{ct, b} = {ft, ft} = {ft, d} = {ft, = 0. (52) 
The first open fermion model is described by the Hamiltonian 

Hit) = H^yg) + H( sys _ env }(t), (53) 
where the system and the interaction are given by 

H(sys) = UJ ftb, 

H [sys _ env) {t) = {(t)$ + £*{t)bc\. (54) 

This model shares all the essential characteristics of a Dirac field interacting with an external 
source which is described by the Lagrangian 



C = ip(h k d k + m)iP + iPZ(t)+Z(t)ip. (55) 
We may find a pair of invariant operators 

/ B(t) \ _ ( u(t) v(t) \fb\f w^(t)c + w«ct \ 

\&(t) J ~ \v*{t) u*(t)J \ft ) + W +) *{t)c + w^>ft J 1 1 ' 



where 



u{t) - + iu )u + (C(t)C(t) + iu> ^)u(t) = 0, (57) 

v(t) - (|| - ^o)v(t) + (C(t)C(t) ~ = 0, (58) 

iw ( -\t) - C(t)u(t) = 0, (59) 

iw {+) (t) + C(t)v(t) = 0. (60) 



Since Eqs. ( pq) and (pOf) are the complex conjugate of Eqs. ( p7|) and (|59f) , respectively, 
solutions can be found of the form 
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v (t) = v*u*{t), w {+) = v*w(~>(t), (61) 
where Vq is a complex number. We further impose the anticommutation relation for all times 

{£(*),#(*)} = 1, (62) 

which leads to the condition 

(u*(t)u(t) + w(->(t)w(-\t)) (l + v*v ) = 1. (63) 

Under the assumption that the external source provides a stochastic force with the correla- 
tions functions 

(C(W*)2 f >(0 = /(*'"*). (W) 

we find the expectation value of the Hamiltonian with respect to the initial thermal equilib- 
rium 



(H(sys))(ther) 



1 V%V 



■(l + 2w { - ) *(t)w { -\t)) 



eP*> + 1 1 - v$v 

+ f dt' f dt"f{t' -t")u*{t')u{t") 



u*(t)u(t)(l-v* Q v ) 

(65) 

The next open fermion model is motivated by a Dirac field interacting with a scalar field 

£ = ip{i^ k d k + rnjif) + (f>(t)tfjijj. (66) 

Again without loss of generality we confine our attention to the single-mode quadratic Hamil- 
tonian 

H{t) = H( sys )(t) + H( S ys-env) (t) , (67) 

where the system and the system-environment Hamiltonians are, respectively, 

H (sys) (t)=rt D \t)(&b-d$), (68) 
H [sys . env) (t) = rt °\t)ft$ - n {0c >(t)M+rt° B \t)Pd-rt 0e >(t)b£. (69) 

We assume = Q( D \ so that the Hamiltonian ( W7h becomes the unitary operator and 

its evolution preserves the unitarity. 



In Ref. |TT| the following type of invariant operators are introduced 



B(t) = u B {t)b + u%{t)ti + Z g(t)d + z^{t)d\ 

bit) = u D {t)h + u + D {t)h ] + z D (t)d + z +(t)d\t). (70) 

In fact, the algebra with the operators in Eq. (|70| ) is closed under the commutation relation 
because of the identity between the commutator and the anticommutators 

[0 X , 6 2 6 z ] = {6 U d 2 }d 3 - 6 2 {d u a 3 }. (7i) 
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The LvN equation leads to the following vector equation 



o 



dt 



u+(f) 
z-{t) 



To simplify Eq. ( [72|) two column vectors are introduced 



LT(t) 



in terms of which it can be rewritten as 



, Z(t) 



.±(U(t)\ /n^ffi Mp(t) \(U(t)\_ 

l dt\Z(t)) + { M G (t) ^ D \t)aJ{z(t)J U - 





fu 


-(f) \ 


n(°-)(t) 




"(*) 





z~ 


"(f) 




w 




fz-{t) + z 






\z-{t)-z 






(U(t)\ 
\Z(t)J 


0. 





(72) 



(73) 



(74) 



Here, 



Mpit) = A<®(t)I + A«(t) ffl + A®(t)a 2 + A( 3 )(t)a 3 , 
M G (t) = -A®(t)I + A^(t)at - A®(t)a 2 - A®(t)a 3 , 

where tXj are Pauli spin matrices and 

A<°)(t) = -^(fi ( ° e) (f) - ^ ( ° e) *(f)), A (1) (f) = ^(^ ( ° e) (f) + ^ (0e) *(f) 



(75) 



(76) 



Further, it was shown in Ref. |TT[ that the anticommutation relations are consistently 
satisfied at all times 



{B(t),B\t)} = {D(t),b\t)} = l, 

{B(t),B(t)} = {B\t),&(t)} = {D(t),D(t)} = {&(t),D\t)} = 0, 



(77) 



Now, B a (t),B^(t) and D a (t), D^ a (t) play the role of the annihilation and creation opera- 
tors for time-dependent fermion system. And the number operators for the particle and 
antiparticle are defined by 



N B (t) = B\t)B{t), 
N D {t) = D\t)D(t), 

whose number states provide the exact quantum states according to Eq. 

N B (t)\n B ,t) = n B \n B ,t), 
N D (t)\n D ,t) = n D \n D ,t), 



(78) 



(79) 



where n B ,nn = 0, 1. 
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V. CONCLUSION 



In this paper we have applied the recently introduced Liouville-von Neumann approach 
to open quantum systems of bosons and fermions. It is based on physically well-grounded 
assumptions. Firstly, the quantum evolution of all systems, either time-independent or time- 
dependent, obeys the Tomonaga-Schwinger equation or the functional Schrodinger equation. 
Secondly, the statistical nature of systems is determined by the density operator. Technically 
the Liouville-von Neumann equation is solved to yield invariant operators which determine 
not only the Hilbert space of exact quantum states but also the density operators. Therefore, 
the Liouville-von Neumann approach describes the nonequilibrium quantum evolution of 
open systems. Moreover, the Liouville-von Neumann approach can also be applied to open 
fermion systems. In summary, the Liouville-von Neumann approach is a unified method that 
can be applied both to the equilibrium quantum evolution of time-independent systems and 
to the nonequilibrium quantum evolution of time-dependent boson and fermion systems. 

The author would like to thank F.C. Khanna, K.-S. Soh and J.H. Yee for collaborations 
and useful discussions. This work was supported by KRF under Grant No. 1998-001-D00354. 
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